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Abstract 



We consider long-baseline neutrino oscillations in the framework of two 
schemes with mixing of four massive neutrinos which can accommodate all 
the existing indications in favour of neutrino mixing. Within these schemes, 
we derive bounds on the oscillation probabilities and the CP-odd neutrino- 
antineutrino asymmetries in long-baseline experiments. Using the results of 
short-baseline neutrino oscillation experiments, we obtain rather strong upper 
bounds on the long-baseline probabilities 1 — Pp^^l and P^ B ^}y Neverthe- 

less, the projected sensitivities of the MINOS and ICARUS experiments are 
better than our bounds. We also show that there are no corresponding con- 
straints fort^— >v\ and'i/^— *v T long-baseline oscillations and that the CP-odd 
asymmetry in the latter channel can reach the maximal value 2/3\/3 allowed 
by the unitarity of the mixing matrix. Some schemes with mixing of three 
neutrinos are also considered. 
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I. INTRODUCTION 



The search for neutrino oscillations remains a central issue of neutrino experiments. In 
several short-baseline (SBL) experiments, which are sensitive to relatively large values of 
the neutrino mass-squared difference 5m 2 (5m 2 > 0.1 eV 2 ), different oscillation channels are 



(-) 



(-) 



(-) 



investigated at present: V^v e (LSND 0, KARMEN @, CCFR @, NOMAD § 
(CHORUS §, NOMAD § and CCFR §). The long-baseline (LBL) neutrino oscillation 
experiments which will operate in the near future will be sensitive to much smaller values of 
Sm 2 (8m 2 > 10~ 3 eV 2 ). Long-baseline experiments with reactor antineutrinos are starting to 
take data (CHOOZ ||) or are under preparation (Palo Verde ||). Long-baseline experiments 
with accelerator neutrinos (KEK-Super-Kamiokande ||10|| , Fermilab-Soudan [11|, CERN- 
Gran Sasso [Q) are under preparation. 

No indications in favour of neutrino oscillations have been found in many experiments 
with terrestrial neutrinos which have been done in the past (see the reviews in Ref. [P~3]]). 
On the other hand, such indications have been found in all solar neutrino experiments 
(Homestake, Kamiokande, GALLEX and SAGE WW). The suppression of the detected 
event rates with respect to those predicted by the Standard Solar Model (SSM) |~5| can 
be explained with neutrino mixing. In the case of resonant MSW transitions ||16|| , it was 
found |17| that the oscillation parameters Sm 2 and sin 2 2t? ($ is the mixing angle) have the 
following values: 



3 x 1(T 6 < 5m 2 < 1.2 x KT 5 eV 2 



4 x 1(T 3 < sin 2 2$ < 1.2 x 10 -2 . 



The solar neutrino puzzle can also be solved by invoking vacuum neutrino oscillations |T8 
in which case Sm 2 ~ 10 _10 eV 2 . 

Indications in favour of z/ M — >u x oscillations (x ^ ^i) have been found also in the 



Kamiokande |T3, 1MB 001 and Soudan ETJ atmospheric neutrino experiments. From the 



analysis of the Kamiokande data the following allowed ranges for the oscillation parameters 



were obtained fi~9| : 



5 x 1(T 3 < 5m 2 < 3 x KT 2 eV 2 
7 x 1(T 3 < 5m 2 < 8 x KT 2 eV 2 



0.7<sin 2 2tf <1 
0.6 < sin 2 2tf < 1 



;i.2) 
;i.3) 



Finally, indications in favour of — > v e oscillations have been found recently in the 
LSND experiment |IJ, in which antineutrinos originating from the decays of /i + 's at rest 
were detected. From the analysis of the data of this experiment and the negative results of 
other SBL experiments (in particular, the Bugey [^] and BNL E776 [23| experiments), it 
follows that 



0.3 < Sm 2 < 2.2 eV 2 



10 -3 < sin 2 2$ < 4 x 10~ 2 . 



;i.4) 



All the above-mentioned hints for neutrino oscillations can be accommodated in schemes 
with four neutrinos. These schemes require a sterile neutrino field (a = s) in addition to the 
usual flavour fields with flavours a = e, fi, r. The specific nature of such 4-neutrino schemes 



has recently been derived from the data f24],|25| . 
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The purpose of the present paper is to consider LBL neutrino oscillations in the light 
of the 4-neutrino schemes favoured by the data. In this context we will discuss the bounds 
on the LBL oscillation probabilities and on the CP-odd neutrino-antineutrino asymmetries 
that can be obtained from the results of the SBL oscillation experiments. We will also make 
an excursion to some 3-neutrino schemes. 

Let us stress that CP violation is one of the important problems which can be tackled 
in LBL experiments. If CP is violated in the lepton sector, then the neutrino mixing matrix 
is complex and the probabilities of u a —>■ Vp and u a — ► Vp transitions (a ^ (3) are different. 
The observation of CP violation in neutrino oscillations will be very important for the 
understanding of the nature of neutrino mixing and the nature of CP violation. In the 
framework of mixing of three massive neutrinos, possible effects of CP violation in LBL 
experiments were considered in Refs. [p6|-|2"8|| . In the present paper we will derive bounds on 



the CP-odd asymmetries in LBL experiments in the framework of the schemes with mixing of 
four massive neutrinos which take into account all the existing neutrino oscillation data. We 
will present general methods for obtaining such bounds from the results of SBL experiments. 
These methods will be applied to some mixing schemes of three neutrinos as well. In the 
context of the 4-neutrino schemes, we will show that sizable CP-odd asymmetries can be 

expected only in i^— >W transitions. 

The plan of the paper is as follows. In Section [TI| we will review the basics of neutrino 
oscillations and CP violation. The 4-neutrino schemes favoured by the data and the relevant 
LBL formulas are introduced in Section |T|. Section |V| is devoted to the discussion of the 
bounds on the LBL oscillation probabilities which follow from the 4-neutrino schemes and 
the SBL data. The same is done for the CP violation parameters in section [V|. Finally, we 
make a digression to 3-neutrino schemes in Section |VT| and we formulate our conclusions in 
Section [VI]] . The three Appendices contain some derivations and discussions of the bounds 
used in the main body of the paper. 



II. NEUTRINO OSCILLATIONS AND CP VIOLATION 

In accordance with the neutrino mixing hypothesis (see, for example, Refs. P^|-|3T|), a 
left-handed neutrino field v a ^ is a mixture of the left-handed components VkL of the (Dirac 
or Major ana) fields of neutrinos with definite masses m k : 

v<xl = ^ Uak VkL with a = e, fJ,, T, s, . . . (2.1) 

k 

where U is the unitary mixing matrix and v s l a sterile neutrino field. The mixing in Eq. ( [2.1[ ) 
implies that the probabilities of transitions of neutrinos (antineutrinos) with momentum p 
at a distance L of the neutrino detector from the neutrino source are given by 

2 

(2-2) 

2 

(2.3) 



>v 



Yl u ? k u *k ex p 



Am 2 kl L 



2p 



ak exp -% 



Am 2 kl L 



2p 
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where Am|j = ml — m\ (we take m\ < rri2 < . . .). From Eqs. (|2.2| ) and ( |2.3| ) it follows that 
the transition probabilities of neutrinos and antineutrinos are connected by the relation 

Pua^Up Pvfj^Va ■ (2-4) 

This relation reflects CPT invariance. 

If CP invariance in the lepton sector holds, then there are phase conventions such that 
in the case of Dirac neutrinos we have 

Uak = Kk , (2-5) 

whereas in the case of Majorana neutrinos we have 

U ak = -U* ak Vk (2.6) 

where r] k = ±z is the CP parity[] of the Majorana neutrino with mass m k (see, for example, 
Ref. |30|). It is obvious that the CP parities rjh do not enter in the expressions for the 
transitions amplitudes. Hence, in both the Dirac and Majorana cases, CP invariance implies 



that [33 



Let us introduce the CP-odd asymmetries 

Da;f3 = Pv a ^v fJ - Pp a ^u ■ (2.8) 

From CPT invariance it follows that 

D a -,p = -D p , a . (2.9) 
Furthermore, from the unitarity of the mixing matrix we have 

£A» ; /j = 0. (2.10) 

We observe that in the case of transitions among three flavour states (u e , z/ M , u T ) the CP 



asymmetries satisfy the relations pi 



D e;fl = D T . e = , (2.11) 

which follow from Eqs. (|2.9| ) and ( 2.10|) . 

In the general case of mixing of an arbitrary number of massive neutrinos, the asymme- 
tries are given by 

D a , p = 4 ^ Im [U aj U; 3 U* k U pk ] sin . (2.12) 

k>j p 

Therefore, CP violation in the lepton sector can be observed only if at least one of the 
oscillating terms in the transition probabilities does not vanish due to the averaging over 
the neutrino energy spectrum and the size of the neutrino source and detector. 



1 The CP parities of Majorana neutrinos could be important for neutrinoless double beta-decay; for 
example, if the v k s have different CP parities, their contributions to the amplitude of neutrinoless 
double-beta decay could cancel each other [32|. 
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III. FOUR MASSIVE NEUTRINOS 



All existing indications in favour of neutrino oscillations can be accommodated by a 
scheme with mixing of four massive neutrinos p3,|36,24|,25| . In Ref . |24] we have shown that 
from the six possible spectral schemes of four massive neutrinos, which correspond to three 
different scales of mass-squared differences Am^, only two schemes are compatible with the 
results of all experiments. In these two schemes the four neutrino masses are divided in two 
pairs of close masses separated by a gap of ~ 1 eV: 



(A) 



atm solar 



LSND 



solar atm 

and (B) m x < m 2 <C m 3 < m 4 . 

s . ' 

LSND 



(3.1) 



In scheme A, Am^ is relevant for the explanation of the atmospheric neutrino anomaly and 
AmJ 3 is relevant for the suppression of solar z/ e 's. In scheme B, the roles of Am^ and Am^ 
are reversed. 

We will consider first SBL neutrino oscillations in the framework of these two schemes. 
Let us note above all that it is impossible to reveal the effects of CP violation in these 
experiments which are sensitive to Am| x > 0.1 eV 2 . Furthermore, both schemes A and B 
give the same oscillation probabilities. Therefore, it is impossible to distinguish A from B 
in SBL oscillation experiments either. In fact, since in these experiments 



Am^ L 



< 1 



the probabilities of u a 



2p 
up and P, 



p(SBL) 

(-) (-) 



and 



Am| 3 L 
2p 



< 1 



(3.2) 



A 



vp transitions are equal and are given by |36 
Am 2 L 



a;/3 



cos ■ 



2p 



(3.3) 



with Am 2 = Am^ 



m\ — m\ and the oscillation amplitude 



A 



a;f3 



k=l,2 



fc=3,4 



(3.4) 



The survival probabilities of neutrinos and antineutrinos are given by 



p(SBL) 

(-) (-) 



(SBL) 

) 



1 



B a - a I 1 



Am 2 L 



cos ■ 



2p 



(3.5) 



where 



fi+a \fc=l,2 / \ fc=l,2 / 

= 4(5^ \u ak \A (1- \ u <* 

\fe=3,4 / V fc=3,4 



(3.6) 
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Long-baseline neutrino oscillation experiments are planned to be sensitive to the "atmo- 
spheric neutrino range" 10 _3 eV 2 < Am 2 ,^ < 10~ 1 eV 2 . In scheme A, the probabilities of 
v a and v a — > vp transitions in LBL experiments are given by 

2 



p(LBL,A) 



-,(LBL,A) 



Upi U* al + U P2 U* 2 exp 



U* pi U al + U* 02 U a2 exp 



(32 ■ 



. Am 2 -! L 
2p 



.A.ml x L 



E 

k=3A 



2p 



y~] Up k u ak 

fe=3,4 



(3.7) 



(3.8) 



These formulas have been obtained from Eqs. fl2.2p and (|2.3|) , respectively, taking into account 
the fact that in LBL experiments Am% 3 L/2p <C 1 and dropping the terms proportional to 
the cosines of phases much larger than 2n (AmL-L/2p 3> 2tt for k = 3, 4 and j = 1, 2). Such 
terms do not contribute to the oscillation probabilities averaged over the neutrino energy 
spectrum. The transition probabilities in scheme B ensue from the expressions ( |3.7| ) and 
fl3.8|) with the change 



1,2^3,4 



(3.9) 



From Eqs.0, Q) and Q it follows that the LBL CP-odd asymmetries £>i^ L) in 
the schemes A and B are given by 



n (LBL,A) 



r(A) • 



2p 



(LBL.B) 

a;/3 



( B) . Am^L 



with the oscillation amplitudes 

= 4 Im [C/ al C/^ EfcJ , if > 



Alm[U a3 U; 3 U* a4 Up,] . 



(3.10) 
(3.11) 

(3.12) 



l af3 — ^ ■ L111 L w « i w a2 , 

We will now discuss the constraints on the neutrino oscillation parameters which follow 
from the existing results of SBL experiments. In the disappearance SBL reactor and acceler- 
ator experiments no indication in favour of neutrino oscillations was found (see the reviews 
in Ref. |R|). At any fixed value of Am 2 in the range 

0.1 eV 2 < Am 2 < 10 3 eV 2 , (3.13) 

from the exclusion plots of the Bugey |E| v e disappearance experiment and of the CDHS 
[37] and CCFR |38j v ^ disappearance experiments we have 

B a ., a <B° a . a (a = e,/i). (3.14) 
Let us define the quantities c a (with a = e, /x, r) in the two schemes A and B as 

(A) c a = \ U ^? , (3.15) 



k=l,2 



(B) 



Cry — 



/] \U a k\ 



(3.16) 



fc=3,4 
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Taking also into account the results of the solar neutrino experiments and the results of the 
atmospheric neutrino experiments, in both schemes A and B, the quantities c e and c M are 
constrained by p4 



where 



The values of a" and 



c, < a c 



c M > 1 - a M , 



B9, 



la 



e,fi) 



(3.17) 



(3.18) 



which have obtained from the exclusion plots of the Bugey |f22 
and CCFR [38|] experiments are given in Fig.l of Ref. []3"9|1 . For values of Am 2 in 



CDHS | 
the range ( |3.13[ ) a° e is small 



£ < 4 x 1(T 2 ), and a° is small for Am 2 > 0.3 eV 2 (oj < KT 1 ). 



In the following we will use also the limits on the amplitudes A a; p which can be obtained 
from exclusion plots of the BNL E734 f4|, BNL E776 H and CCFR 



v.. 



(-) 

*v e appearance 



experiments and of the FNAL E531 g| and CCFR 



(-) 



(-) 



appearance experiments: 



.4 



a;/3 



< A 



(3.19) 



with a = fi and (3 — e, r. The results of the LSND [|TJ — > v e appearance experiment will 
be also taken into account (in particular, the allowed range ( |1.4| ) of Am 2 ). 

Since scheme B emerges from scheme A by the substitution ( |3.9| ) and since we will derive 
bounds on the LBL oscillation probabilities Prf'^w' , P ( _^ B ^f and on the CP-odd parameters 

Ia3i lap as functions of A a -,p, c a and c^, it is evident that such bounds apply equally to 
both schemes A and B by virtue of the definitions ( |3.4j ), ( |3.15| ) and ( |3.16| ). Consequently, 
when dealing with such bounds we will omit the superscripts A, B indicating the specific 
scheme. 



IV. CONSTRAINTS ON THE LONG-BASELINE PROBABILITIES 

We will consider first the limits on the LBL oscillation probabilities which can be obtained 
from the results of the SBL oscillation experiments. The Cauchy-Schwarz inequality implies 
for scheme A that 



Yl U P k U *k ex P 



fc=l,2 



Am 2 kl L 



2p 



< C a C{3 



(4.1) 



Using this inequality and the definition ( |3.15| ) of c a , we find from the LBL probabilities in 
Eqs.( |3.7| ) and ( |3.8| ) that the survival probabilities Pi^^ Q and Pf£^f£ a are bounded by 



(i - c a f < ifff* < c 2 + (i - Ca y . 



-) (-) 



(4.2) 



As explained at the end of the last section, these bounds are scheme-independent. In order to 
obtain bounds on the LBL transition probabilities pj^tj^ and Ppl^Jg, we take into account 
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the definition ( |3.4j ) of A a] p and the inequality (}4.1|). When inserted into Eqs. (|3.7|) and 
they imply 

- A a , p < i^O < c a c fi + - A a , p . (4.3) 

The bounds Q4.2Q and ( |4.3| ) are the basis of the following considerations for the different 
oscillation channels in LBL experiments. 

The smallness of c e (see Eq.( [3.17| )) implies that the electron neutrino has a small mix- 



ing with the neutrinos whose mass-squared difference is responsible for the oscillations of 
atmospheric and LBL neutrinos (z/i, v-i in scheme A and z/3, v± in scheme B). Hence, the 
probability of transitions of atmospheric and LBL electron neutrinos into other states is 
suppressed. Indeed, taking into account the constraint c e < a°, the lower bound on Pj£^p e 
and the upper bounds on Pp^} e which we will derive are rather strict. 

Let us discuss first the bounds on the LBL survival probability Pj£^p B . We will compare 
these bounds with the expected sensitivity of the CHOOZ || and Palo Verde |J LBL reactor 
experiments. Taking into account the constraint ( |3.17| ) on c e , Eq.( |4.2| ) implies that in both 
schemes A and B 

1 - J^j < al (2 - a°) . (4.4) 
The curve corresponding to this limit obtained from the 90% CL exclusion plot of the Bugey 



22] experiment is shown in Fig.|l] (solid line). The expected sensitivities of the LBL reactor 
neutrino experiments CHOOZ and Palo Verde are also shown in Fig.|l] as the dash-dotted 
and dash-dot-dotted vertical lines, respectively. These expected sensitivities with respect 
to 1 — -Pj^y have been extracted by us from the figures presented in Refs. ||[| showing 
the sensitivity of the respective experiments in the sin 2 2'd-Sm 2 plane, using the fact that 
for high values of 6m 2 each experiment is sensitive only to the averaged survival probability 
Pp^&i — 1 — \ sin 2 2$. Thus, the vertical lines in Fig.[l| correspond to ~ sin 2 2t? at high 8m 2 . 

Figure |l| shows that, in the framework of the two schemes (|3.1| ) with four neutrinos, which 
allow to accommodate all the indications in favour of neutrino oscillations, the existing data 
put rather strong limitations on the probability of LBL transitions of v e into other states 
(for Am 2 > 3 eV 2 the upper bound for 1 — Pj^y is close to the border of the region of 
sensitivity of the CHOOZ experiment, whereas for Am 2 < 3eV 2 it is much smaller). 

The shadowed region in Fig.|] corresponds to the range (|1.4|) of Am 2 allowed at 90% 
CL by the results of the LSND and all the other SBL experiments. It can be seen that the 
LSND signal indicates an upper bound for 1 — P^^p e of about 5 x 10 -2 , smaller than the 
expected sensitivities of the CHOOZ and Palo Verde experiments. 

Let us stress that, in the framework of the schemes under consideration, the smallness 
of c e is a consequence of the solar neutrino problem. Consider for example scheme A. The 
probability of solar neutrinos to survive is given by 

k=l,2 \ fe=l,2 / 
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where Pj%-} Ve is the survival probability due to the mixing of v e with u% and ^4, depending 
on the small mass-squared difference Am^. From the results of SBL reactor experiments 
it follows that the quantity c e = \U e k\ 2 can be small or large (close to one). In order 

fe=l,2 

to have the energy dependence of the survival probability Pu^ve an d the suppression of 
the flux of solar z/ e 's that are required for the explanation of the data of solar neutrino 
experiments, we must choose a small value of c e . In this case, the survival probability of 
z/ e 's in LBL reactor experiments is close to one. 

Let us now discuss the bounds onV^— >u e transitions in LBL accelerator experiments. We 
will compare these bounds with the expected sensitivities of the KEK-Super-Kamiokande 
(KEK-SK) 0, Fermilab-Soudan (MINOS) and CERN-Gran Sasso (ICARUS) | 



experiments. 

Taking into account the constraints (|3.17|) on c e and (|3.19|) on A^. e , Eq. ( |4.3|) implies that 
in both schemes A and B 

The curve corresponding to this limit obtained from the 90% CL exclusion plots of the Bugey 
Z^j experiment for a° e and of the BNL E734 BNL E776 |23| and CCFR |7[ experiments 



for A°^. e is shown in Fig.^j (long-dashed line). For a comparison, we have shown the expected 
sensitivities of the LBL accelerator neutrino experiments KEK-SK [10|, MINOS and 



ICARUS |12| (the dotted, dash-dotted and dash-dot-dotted vertical lines, respectively). 
These sensitivities have been obtained from the figures presented in Refs. [[10HT2| showing 



the sensitivities of the respective experiments in the sin 2 2'dSm 2 plane with the method 
explained in the context of LBL reactor experiments. 

The conservation of probability and Eq. (f4.2|) lead to a further upper bound: 



P<£t% < 1 - PSS a <c a (2- c a ) (a ?0). (4.7) 



CPT invariance (see Eq. Q2.4j) ) and the fact that the bound (|4.2|) is valid for antineutrinos 
as well give the same upper bound as Eq.( |4.7|) for the opposite transition up — >• v a : 

P^l = P&1% < 1 - P$SL <c a (2- ca) (a?0). (4.8) 

Finally, these two equations hold evidently also for antineutrinos. The solid curve in FigfJ 
represents the limit 

^ < a\ (2 - al) (4.9) 



obtained from Eq. ( |4.8| ) and the constraint ( 3.17| ) on c e . This bound is better than the bound 



1|) for the SBL parameter Am 2 < 0.4 eV 2 . 
The darkly shadowed area in Fig.|2] represents the region allowed by the results of the 
LSND experiment, taking into account the results of all the other SBL experiments. The 
lower bound on i^ 5 ^ , i.e. the left edge of the shadowed region, is determined by the results 
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of the LSND experiment, whereas the upper bound is given by the limits obtained before, 
i.e. by the most stringent of the inequalities ( |4.9|) and (|4.6p represented, respectively, by the 
solid and long-dashed curves. Thus, for each fixed value of Am 2 we have 

where A™™ is the minimal value of A^ e allowed at 90% CL by the LSND experiment. 

Figure ^] shows that, in the framework of the schemes under consideration, the sensitivity 
of the KEK-SK experiment may be not sufficient to reveal LBL z/ M — > v e oscillations, whereas 
the sensitivities of the MINOS and ICARUS experiments are considerably better than the 
upper bound for pj)^} e - ft is interesting to observe that there is also a lower bound on 
this transition probability that follows from the LSND results (see Eq. (|4.10 )). However, this 



lower bound is valid only in the case of LBL neutrino oscillations in vacuum. The corrections 
due to the matter effects in LBL experiments make it disappear (see Section[VTI|) . 

The solid curve in Fig.|2] is at the same time an upper bound on P^fX . This is evident 

from Eq.( [4.7|) . On the other hand, the probability of Uf l —*v T transitions is not constrained 
by the results of SBL experiments. 

Finally, a further upper bound on P^^\ for a ^ (3 is gained from Eq.(O). Since 

A a; p < 4(1 — c a )(l — Cp), we have 

i^fS < c Q cp + (1 - c Q )(l - cp) (a^P). (4.11) 

Obviously, if c a = Cp = or 1 is in the allowed range of these quantities, then this upper 
bound is 1 and thus is trivial. This leaves only a = fi and (3 = e as a non-trivial case, with 

p ( TcU a l+<- 2a K- ( 4 - 12 ) 

The short-dashed curve in Fig.|2| shows this limit with a° e and a° obtained from the 90% CL 
exclusion plots of the Bugey |22|] v e — > u e experiment and of the CDHS [^7] and CCFR [38 



experiments, respectively. For a" C < 1 this bound is about half of that given 
by Eq. (|4.9|) . However, since a° is only small in the same range of Am 2 where A°^. e is small, 
numerically the bound ( 4.12 ) turns out to be worse than the bound ( |4.6| ) (the long-dashed 
curve in Fig.0). 

V. CP VIOLATION IN THE SCHEMES WITH FOUR NEUTRINOS 

As shown in Appendix [FJ, the unitarity of the mixing matrix implies the "unitarity 
bound" 

1^1 < f(c a ,Cp) (5.1) 

where f(x, y) is the continuous function 
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fi = xy for 2(1 — x){l — y) > xy 

f(x,y) = { ' (5.2) 
h = 2[(x + y-l)(l-x)(l- y)} 1 ' 2 for 2(1 - x )(l-y)< xy 

defined on the unit square < x < 1, < y < 1. In Fig.|3| we have drawn a contour plot 
of the function f(x,y), which is helpful for the determination of the maximal allowed value 
for f{c a ,c@) when c a and/or eg are bounded. The dotted line in Fig|| is the borderline 
g(x) = 2(1 — x)/(2 — x) between the regions where / = f\ and / = / 2 . Note that / is 
continuous along this borderline. 

In order to determine the maxima of f(x,y), the following considerations are useful 
(for the details consult Appendix |Cj). Increasing x at fixed y, the function / increases 
monotonously from / = at x = 0, until the straight line yi(x) = 2 — 2x (1/2 < x < 
1) depicted in Fig|3] is reached. There, the value of / is given by / = yy/1 — y. After 
this intersection, the function / decreases monotonously to / = at x = 1. From the 
symmetry f(x,y) = f(y,x), it follows that for fixed x and increasing y the function / 
increases monotonously from / = 0aty = 0to/ = xy/1 — x when the straight line 
2/2 (#) = 1 — x/2 (0 < x < 1) depicted in Fig.||is crossed. After this intersection, / decreases 
monotonously to / = at y — 1. 

The absolute maximum of the function / (see Appendix |C|) lies at the intersection of the 
lines y\ and y% and is given by / max = 2/3^/3 ~ 0.385. Therefore, from the unitarity of the 
mixing matrix we have an absolute maximum for \I a p\: 

%p\ < ^ « 0.385 . (5.3) 



With the help of FigfJ, one can see that Eq.( |5.1|) with the constraints (|3.171) on c e and 
c M implies that 



, / 2 (a°,y 2 (a°)) = a° e (l - a°) i/z for a° > a°/2 , 

VI < { , ( 5 - 4 ) 

fM 1 - al) = 2 [(a° - (l - a°) aj] 1/2 for a° < a°/2 . 

The solid curve in Fig.^ shows the limit \I efl \ < a° ^Jl — a® with a° e obtained from the 90% 



CL exclusion plot of the Bugey [22] v e v e experiment. The dotted curve in Fig .H represents 



the improvement reached with the lower part of Eq.fl5.4p at the values of Am 2 for which 



< a° e /2, with < obtained from the 90% CL exclusion plots of the CDHS P] and CCFR 



38] { v il —*v iJi experiments. 



The bound represented by the solid curve in Fig.|] is valid also for |/ er |, because there is 
no experimental information on c T . 

For \Ifj, T \, again by inspection of Fig.|^, one can see that Eq. (|5.1| ) with the constraints 



( p. 17] ) on c M implies that 

I VI < / 2 (1 - oj, y 2 (l - al)) = (1 - al) yfa . (5.5) 

The solid curve in Fig|| represents the corresponding bound obtained from the 90% CL 
exclusion curves of the CDHS |37j and CCFR |38[] ^v^—^v^ experiments. For Am 2 < 0.3 eV 2 
there are no experimental data and therefore | Wax ~ 0.385 by virtue of Eq.(p.3|). 
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Taking into account the expression (|3.4j) for A a -p, in both schemes A and B we have (for 
the proof of this inequality, see Appendix |A|) 



Iaf3 1 < ^ V Aa;f3 (4 C Q Cfl - A a -0 ) 



Taking into account Eq. fl3.19|) , we obtain 



l a/3\ 



< 



C Q C/J 



for 
for 



A° a .p<2c a c p 
A° a . p >2c a c p 



For |/ eM |, with the constraints ( 3.17| ), the inequality ( |5.7| ) becomes 



(5.6) 



(5.7) 



< 



for 
for 



^ ;e >2a°. 



(5.8) 



The dash-dotted curve in Fig.[| shows the limit (|5.8|) obtained using the 90% exclusion plots 
of the Bugey |22 v e — > v e experiment for the determination of a° e and the BNL E734 ||40|| , 
BNL E776 p3[ and CCFR B v„ — > i/ e experiments for the determination of A . 



Since the constraints (|3.17|) do not put an upper bound on the possible values of c M and 
c T , in the case of |J /JT | the inequality ( |5.7| ) becomes 



(5.9) 



The dotted curve in Fig.|| shows the limit (|5.9| ) obtained using the 90% exclusion plot of the 

v T experiment for the determination of A^. T . 



FNAL E531 ET| and CCFR 



v.. 



The shadowed regions in Figs.|] and |^ correspond to the range ( |1 .4j) of Am 2 allowed at 
90% CL by the results of the LSND and all the other SBL experiments. From Fig.|5] it can 
be seen that, taking into account the LSND signal, |J^ r | could be close to the maximal value 
2/3 \J% allowed by the unitarity of the mixing matrix. 



VI. THREE MASSIVE NEUTRINOS 

It is worthwhile to have a look at LBL neutrino oscillation experiments neglecting some 
of the present hints for neutrino oscillations. It is possible that not all these hints will be 
substantiated in the course of time and it is useful to check which features are actually 
dependent on or independent from them. 

In this Section we consider the minimal scenario of mixing of three neutrinos. We 
will assume that of the two differences of neutrino mass-squared one is relevant for SBL 
oscillations and the other one for LBL oscillations (see also Refs. ||2"6j-|2"8"|). Hence, in this 
section we adopt the point of view that not neutrino mixing but other reasons could explain 
the solar neutrino data. With these assumptions, there are two possible three-neutrino mass 
spectra: 
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LBL 

(I) mi < m 2 < m 3 

S v ' 

SBL 



LBL 

and (II) mi <C m 2 < rm 
s . ' 

SBL 



(6.1) 



In both schemes I and II, Am| x is assumed to be relevant for neutrino oscillations in SBL 
experiments. In this case, the SBL oscillation probabilities depend on |f/ e3 | 2 and {U^l 2 in 



the scheme I [|39| and on \U e \\ 2 and iC^i) 2 in the scheme II ||36|| . There are three regions of 



these quantities which are allowed by the results of disappearance experiments (see Refs. 



(1) 
(2) 
(3) 



\U ek \ 2 > 1 



\U, 



ek 



\u, 



ek 



<a 



\u, 



iik I 2 < ( 'l > 

lik] 2 < "I , 



(6.2) 



lik | 



> 1 



with k = 3 for the scheme I and k = 1 for the scheme IlQ (for the definition of a° e and a°, 
see Eq.flngD). 

The neutrino and antineutrino LBL oscillation probabilities in scheme I are given by 



p(LBL,I) 



,(LBL,I) 



Ufsi U* al + U P2 U* 2 exp 



Upx U a i + U*p 2 U a2 exp 



2p 

. Am?,! L 
2p 



Up 3 \ 2 \U, 



a3\ 



2 \TT I 2 
\ u a3\ 



(6.3) 
(6.4) 



with Eqs.(p.7|) and (p.8|), it is obvious that the 



From the comparison of Eqs. ( |6.3| ) and 

CP-odd asymmetries D^ h,V> are given by the same formula ( |3.10| ) as in the 4-neutrino case 



(With ^ *a0 

expressions ( |6.3D and (|6.4j ) with the cyclic permutation of the indices 



1^1). The transition probabilities in the scheme II can be obtained from the 



1, 2, 3->2, 3, 1 . 



(6.5) 



Therefore, as in the case of the schemes A and B for four neutrinos, the bounds on the LBL 
oscillation probabilities and the CP-odd asymmetries are the same in the three neutrino 
schemes I and II. 

The bounds on the LBL oscillation probabilities E^ B f\ and on the CP-odd parameters 

I a/ 3 derived in the Appendices for the four-neutrino schemes ( |3.1| ) are valid also in the case 
of mixing of three neutrinos: the demonstrations in the four-neutrino case A (B) can be 
applied to the three neutrino case I (II) if we put U a 4 = (U a i = and change the indices 
2, 3, 4 — > 1, 2, 3) for all a = e, fi, r. 



For a comparison, the schemes I, II and the regions 1, 2, 3 are called hierarchies II, I and regions 



A, B, C, respectively, in Ref. [26| 
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It is obvious that, with A a -^ = 4|[/g3| 2 |[y Q 3| 2 , the same bounds on i?^ 3 ^ arise for 

Va—*Up 

a = (3 and a ^ (3 as given by Eqs.([4.2[) and ([4.3|). Since in the 3-neutrino case the CP-odd 



asymmetries in different oscillation channels are connected by Eq.( 2.11 ), we have 



lea = Lt = Ire • (6.6) 



A few remarks on the unitarity bound ( |5.1|) are in order. It is true that from the unitarity of 
the 3x3 mixing matrix we have A a .p = 4(1 — c a )(l — c^) and c a + cp > 1, but, nevertheless, 
the distinction defined by Eq. (|B6|) has to be maintained. Therefore, also the unitarity bound 
is upheld with the addendum that c a and cp can only vary within the section of the unit 
square defined by c a + > 1. Since the point c a = c@ = 2/3 fulfills this condition, the 
absolute maximum ll^max — 2/3y / 3 of the 4-neutrino case extends its validity to three 
neutrinos^. 

In the following we will give LBL bounds for each of the regions ( |6.2|) , along the lines of 
the previous 4-neutrino sections. 

Region 1. With respect to SBL and LBL neutrino oscillations, the 3-neutrino schemes 
I and II in Region 1 correspond to the 4-neutrino schemes A and B, respectively, with the 
same bounds on P^ B ^ (Eq.flj) and Fig.0), i^ B £> (Eqs.Q, Q, ( ^T2|) and Fig.@) and 



I ejJ \ (Eqs. (|5T4T ), (|5\8| ) and Fig.f|). Since \I^ T \ = |I e J, the stringent bounds on |J e J given in 



Fig.f| are valid also for \I^ T \. 

For completeness, we want to mention that there is a change in the upper bound for 
P ( l^f \ in going from four to three neutrinos: taking into account the inequality c e + c„ > 1, 



(-) (-) 



we have c M > 1 — min(a°, a°) and Eq. (|4.12|) improves to 

P? B fi < a° e + (1 - 2a° e ) min(a°, aj) . (6.7) 



For a e < a" this bound is slightly more stringent than that given by Eq. (|4.9| ), but the 
improvement is negligible for a® <C 1. 

Region 2. Actually, this Region is excluded by the results of the LSND experiment (see 
Refs. pP|J3B|f2^] ] ) . The reason is that (in combination with other data) the upper bound 

A Ke <4a° e al (6.8) 

is too restrictive to be compatible with the LSND data. In spite of this evidence, let us 
discuss the bounds on the LBL probabilities in this Region. 

The restrictions c e > 1 — a°, c M > 1 — a° and the unitarity of the mixing matrix imply 

V < a° e + a°. 



that c r is small: c r = 2 — c e — c M < a° e + a°. From Eq. (|4.3|) it follows that the probabilities 
of { v il — >v T andVg— >v T transitions in LBL experiments are confined in the range 



3 This value is 4 times the maximal value of the Jarlskog parameter J [42| for CP violation in the 
Kobayashi— M askawa matrix, | «/| max = l/6\/3. 
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1 • - ^(LBL) . 1 , , _0 , _0 



7 A aiT < < - A a;T + + a£ (a = e, //) , (6.9) 



whereas for the probability of v^— >v e transitions we have only the lower bound 

\A Ke <P^ B V. (6.10) 
The inequality (|4.8|) yields the additional upper bounds 

^T( L) , < K + <) t 2 - «° - aj) (a = e, . (6.11) 

The bounds on CP violation can be derived with the methods described in the Appen- 
dices. They are given by the oscillation amplitude bound 



\U<l^e (4-^ ;e ) (6-12) 

and the unitarity bound 

|/e„l < / 2 (1 - a°, 1 - oj) = 2 (l- a O- a o) . (6.13) 

Both are of similar order of magnitude and less restrictive than the bounds in Region 1. 
Since \I^ T \ = \I e ^\, these bounds are valid also for \I^ T \. 

Region 3. In this Region, where c e > 1 — a° and c M < a°, the atmospheric neutrino data 
cannot be explained in the framework discussed here. The reason is that 

> (1 - a ,)' (6.14) 

and this is incompatible ]23j with the atmospheric neutrino anomaly. The LBL transition 
probabilities of muon neutrinos are confined by 

7 A^p < < - A, ]P + aj ((3 = e, r) , (6.15) 

whereas for ^v e — >v T transitions there is only the lower bound 

-A e , T <P^. (6.16) 

The inequality ( |4.7|) , which is a consequence of probability conservation, leads to 

P^} } < aj (2 - aj) 09 = e,T). (6.17) 

Furthermore, taking into account the inequality c e + c M > 1, we have c e > 1 — min(ag,a°) 
and Eq.( (4.12|) improves to 
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iffS <<+(!- 2aJ) min(a°, a J) . (6.18) 



For a° -C a° 1 this bound is about half of that given by Eq. fl6.17p . 
Finally, our methods for obtaining bounds on CP violation yield 



\I e »\<-y/A% e (±a°-A% e ) (6.19) 



and 




for a° > a°/2 , 
i£K] 1/a for a°<a°/2. 



(6.20) 



The oscillation amplitude bound is more stringent in this case. From Eq. (|6.6|) it follows that 
the bounds ( 6.19|) and ( |6.20 ) are valid also for the parameter \I^ T \ that characterizes the 



CP-odd asymmetry in the — > v T channel. 

The differences in the bounds on the LBL probabilities are marked and could thus serve 
to distinguish between the three different Regions in the three neutrino case. Of course, in 
the experiments discussed here the four neutrino case (schemes A and B) is indistinguishable 
from the three neutrino case with Region 1. The above-mentioned distinctions could also 
serve as a cross-check for present hints of neutrino oscillations. 



VII. CONCLUSIONS 

At present there are three experimental indications in favour of neutrino oscillations 
which correspond to three different scales of neutrino mass-squared differences: the solar 
neutrino deficit, the atmospheric neutrino anomaly and the result of the LSND experiment. 
These indications and the negative results of numerous short-baseline neutrino experiments 



can be accommodated in two schemes (A and B) with mixing of four massive neutrinos J24 
In this paper we have presented a detailed study of the predictions of the schemes A and 
B for long-baseline experiments. Using only the results of existing experiments, we have 
obtained bounds on probabilities of different transitions in long-baseline experiments. 

The schemes A and B give completely different predictions for neutrinoless double beta 
decay and for neutrino mass effects in experiments for neutrino mass measurements by the 



tritium method [23]. They lead, however, to the same bounds on long-baseline oscillation 



probabilities. In addition, all the bounds that we have derived apply for neutrinos as well 
as antineutrinos. 

We have shown that the results of the short-baseline reactor experiment put rather 
severe bounds on the probability 1 — P D \^ e of v e transitions into all possible other states 
in long-baseline experiments. If the Am 2 relevant in short-baseline oscillations is bigger 
than about 3eV 2 , the bound on 1 — Pp^& e is slightly higher than the sensitivity of the 
CHOOZ experiment, allowing some possibility to reveal neutrino oscillations in this channel. 
However, the results of the LSND experiment favour the range 0.3 < Am 2 < 2.2 eV 2 . We 
have shown that in this range the upper bound for the quantity 1 — Pp^) lies between 10~ 2 
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and 5 x 10 2 (see Fig.|T]) and thus below the sensitivity of CHOOZ and Palo Verde. On the 
other hand, there is no restriction on i^ BI ^ 

The probability Pr^YK of z/ M — >v e transitions is severely constrained in the schemes A 

and B by the results of short-baseline reactor and accelerator experiments (see FigfJ). The 
sensitivity of MINOS and ICARUS is well below the upper bound for Pzff\, whereas the 

sensitivity of the KEK-SK experiment seems to be insufficient. There is also an upper bound 
on long-baseline { v e —±v T oscillations (solid curve in Fig.^), which is nearly as tight as the one 

for'z/^— >u e transitions. On the other hand, the long-baseline ( z/^—>z/| channel is unconstrained. 

We have obtained bounds on LBL transition probabilities in the case of the neutrino 
mass spectra Q3.1|) , which are implied by the results of the solar, atmospheric and LSND 
experiments. If the LSND data are not confirmed by future experiments, but nevertheless 
there is a mass (or masses) approximately equal to 1 eV providing an explanation for the 
hot dark matter problem, then the neutrino mass spectrum can be different from the spectra 
A and B in Eq.( |3.1| ). The natural neutrino mass spectrum in this case is hierarchical and in 
this case the bounds that we have obtained in this paper are not valid. 

In the framework of neutrino mixing in schemes A and B, we have also derived constraints 
on the parameters I a p that characterize the CP-odd neutrino-antineutrino asymmetries 
in long-baseline experiments. We have developed methods for deriving upper bounds on 
the parameters I a p from the data of short-baseline experiments, which can be applied to 
different schemes. We have shown that CP violation in the — > v e channel is bounded by 
\I efl \ < 10~ 2 . A similar suppression of CP-odd effect takes place in v e — > v T long-baseline 
neutrino oscillations. On the other hand, sizable CP violation can be expected in z/ M — > v T 
oscillations. The CP-odd asymmetry in this channel could be close to its maximally allowed 
value |/ MT |max ~ 0.385, resulting from unitarity of the mixing matrix. 

Summarizing, we would like to emphasize that the results of all neutrino oscillation 
experiments lead to severe constraints on the probabilities for the probabilities for v e dis- 
appearance and ( ^— >v e and ( v e — >v T appearance in long-baseline experiments. Nevertheless, 
even these channels are within the planned sensitivity of the MINOS and ICARUS experi- 
ments. The channels^— >u T and i^— are not constrained at all. Therefore, from the point 
of view of the present investigation, long-baseline muon neutrino beams provide promising 
facilities for the observation of neutrino oscillations. 

Let us stress that we have considered here vacuum LBL oscillations. In conclusion, we 
want to present some remarks about the effects of matter. In the 4-neutrino schemes under 
consideration active neutrinos can transform into sterile states. Therefore, not only charged 
current but also neutral current interactions must be taken into account [16| and the effective 
Hamiltonian in the flavour representation is given by 

H cS = (uM 2 tf + diag {a CC , 0, 0, a NC )) (7.1) 

with 

acc = 2V2G F N eP ^8xlO- 5 eV 2 (—^)(-?—) , (7.2) 

\lgcm _d / VIGeV/ 
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a NC = V2G F N n p ~ ^a C c 



(7.3) 



where M 2 is the diagonal matrix of the squared masses, Gp is the Fermi constant, N e and N n 
are the electron and neutron number density^, respectively, and p is the density of matter, 
which in the Earth's crust is of the order of 3 g cm -3 . The parameters acc and a^c are small 
with respect to the Sm 2 relevant for LBL oscillations. Nevertheless, their contributions for 
some transitions can be important (see |28) for the case of three neutrinos). 

In order to take into account the matter effects in the oscillation probabilities we must 
diagonalize H e f^\. In the case of the schemes (|3.1|) under consideration, the relation acc <^ 



(7.4) 



Am ~ leV is valid and thus, up to terms of the order acc J 'Am , we have 



e 

2p 



where e = diag (ei, . . . , 64), U' — UR and R has the block structure 

r \ R ' 

- 1 R" 



(7.5) 



The matrices R' and R" are 2x2 unitary matrices. Now instead of Eq.(3.7) we obtain for 
the LBL oscillation probabilities 



p(LBL) 



2 

E W ex p (-*§ L ) + £ exp L ) 



(7.6) 



From the Cauchy-Schwarz inequality it is easy to see that in both schemes A and B the 
upper bounds (|4.2|), (|4.11|) and therefore ( |4.12| ) (represented by the short-dashed curve in 
Fig.||) are valid in the case of matter apart from corrections of the order of ac*c/Am 2 . Other 
bounds could be modified by the matter effect. 

However, the upper bound ( |4.4| ) applied to reactor LBL experiments is practically not 
modified by matter effects. In fact, the relevant parameter e± — €3 ~ accL/2p w 0.5 x 
10~ 6 (L/lm), for p = 3 gem -3 , is around 10 -4 for CHOOZ and Palo Verde, whereas in the 
case of MINOS and ICARUS it is not small (a cc L/2p w 0.4). For the KEK-SK experiment, 
matter corrections are still modest because the baseline is roughly a factor of three smaller 
than the baseline of MINOS and ICARUS. Due to matter effects the bounds ( f4.6|) and ( |4.9|) 
on the Up — > v e transition probability depicted in Fig.^ become less restrictive, in particular 
for the MINOS and ICARUS experiments. Note that also because matter induces corrections 
of the order of 10 -2 the lower bound obtained in the vacuum case from the results of the 
LSND experiment (the left borderline of the darkly shadowed area in Fig.^|) disappears and 
the lightly shadowed region in Fig.0 is allowed by the LSND experiment in the case of LBL 



4 iV e ~ N n ~ ^igi w h ere -^A is the Avogadro number. 

5 Here we consider the simplifying approximation of constant N e and ^V n . 
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oscillations in matter. A detailed discussion of matter effects in the case of mixing of four 
massive neutrinos will be published elsewhere [Q . 

For antineutrinos the relevant effective Hamiltonian is obtained from Eq. ( |7. 1\ ) by taking 
U* instead of U and putting a minus in front of acc and a^c- Therefore, matter effects 
can mimic CP-odd asymmetries and at least at the order of 10 -2 both effects are entangled. 
In the framework of the schemes with mixing of three neutrinos the effects of matter in 
long-baseline experiments were considered recently in Refs. pB|-pS|. It was shown that the 
contribution of matter to the CP-odd asymmetries is small for all transitions, except for 
v ii — > v e transitions in the Region 2 of the scheme I which is not compatible with the 



LSND result 39 
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APPENDIX A: DERIVATION OF THE OSCILLATION AMPLITUDE BOUND 



In this Appendix we discuss the derivation of the "oscillation amplitude bound". The 
starting point is the quantity 



I a/3 = 4 Im [U al U* pi U* 2 Up] (a^P), 

which determines the CP-odd asymmetry in the case of four massive neutrinos. 
It is obvious that I a p is invariant under the phase transformation 



(Al) 



77 . _> e l 7j TJ . 



U 



Hi 



e^U 



(A2) 



where the 7^ are arbitrary phases. Thus the elements U a j can be taken to be real. Taking 
into account the definitions ( |3.15| ) and ( p . 1 6[ ) valid in the schemes A and B of Eq.( |3.1| ), 
respectively, we can write 



U 



'3 



with 



j = 1, 2 in the scheme A, 
j = 3, 4 in the scheme B, 



and the orthonormal basis 

e«(fl) = (cos#, sin 1 



sin 9, cos ( 



(A3) 



(A4) 



We expand Upj (with j = 1, 2 in the scheme A and j = 3, 4 in the scheme B) with respect 
to this basis as 



(p) 



(A5) 



p=l,2 



where pi and j»2 are complex coefficients such that 

E = 1 ■ 

p=l,2 



(A6) 



With the help of Eqs.Q-Q we easily find 

I a/3 = 2c a cp sin 26 1 Im(p*p 2 ) = 2 c a |pi | yl — \pi\ 2 sin 2# sin 5 , (A7) 

where 5 is the phase of p*P2- 

The parameter \pi\ is connected to the oscillation amplitude A a - t p and the parameters 
c a , eg. In fact, from Eqs. ( |A3|) and (|A5| ) we have 





4 


E u °* u h 






3=1,2 




4 


E ^ 






j=3,4 



4c Q C/3 |pi| 2 in the scheme A 



4 c Q C/3 1 pi I in the scheme B . 



(A8) 
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Hence, in both schemes A and B, we have \pi\ = -v/ A a .p/ 'Ac a cp. Inserting this value in 
Eq.(|A7|), we obtain 



I a/3 = 2 y A a;/3 (4 c a cp - A a .p) sin 29 sin 5 
and thus we arrive at the "oscillation amplitude bound" 



1 1 a/3 1 < £ \l A a .p (4 C a Cp - A a .p i 



(A9) 



(A10) 



Let us stress that this derivation is based only on the obvious inequality 

| sin 29 sin S\ < 1 . (All) 
Since c a , c@ and A a .p do not restrict 9 and 5, the bound ( |A10| ) is the optimal one. 

APPENDIX B: DERIVATION OF THE UNITARITY BOUND 



Up to now we did not use the unitarity of the mixing matrix. Taking this fact into 
account will allow us to obtain an upper bound on \I a p\ depending solely on c a and Cp. 
The unitarity of the mixing matrix tells us that 



(Bl) 



E u *i u h = - E u ^ u h ■ 

J=l,2 j=3A 

This relation allows to write the oscillation amplitude A a -p in the two forms of Eq.( |3.4| ). 
Using the Cauchy-Schwarz inequality, one can see that 



Aa;(3 = 4 



E U <*3 U h 
i=3,4 



< 



in the scheme A, and 



Aa;(3 = 4 



E u °i u h 

i=i,2 



< 



'(em)(e 

\i=3,4 / \i=3,4 

'(em)(e 

\i=l,2 / V=l,2 



4(l-c Q )(l- C/3 ) 



vB2) 



u Pj \ 2 



4(l-c Q )(l- C/3 ) , (B3) 



in the scheme B. Hence, in both schemes A and B the oscillation amplitude A a .p is bounded 
by 



A a ,p < 4 (1 - Ca) (1 - c p ) 



(B4) 



The right-hand side of the inequality (|A10|) , as a function of reaches its maximum, 

Cq,C/3 5 at 



(A a -,p)o = 2c a cp . 



(B5) 



Consequently, if the condition 
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2 (1 - c Q ) (1 - C/3 ) >c a c p (B6) 

is satisfied, the upper bound (]B3 ) on A Q;/ g is larger then (y4 Q;/ g) . In this case we have 

1 4/3 1 < c a c p . (B7) 

If the condition ( |BT)| ) is not fulfilled, the upper bound is smaller than (A a;/3 ) and has 
to be inserted for A a .p into Eg. ( |A10| ), leading to 



\I«p\ < 2 y (Ca + C/3 - 1) (1 - Ca) (1 - C/9 ) . (B8) 

Thus, we arrive at the "unitarity bound" 

\I a p\<f{c a ,c p ), (B9) 

with the function 

f fx=xy for 2(1 - y) > , 

f(x,y) = { (BIO) 
[ /a = 2[(x + y - 1)(1 - x)(l - y)] 1 / 2 for 2(1 -x)(l-y)<xy, 

defined on the unit square 0<x<l,0<y<l. The function 

, , 2(1 -x) . 
M = -^f- (Bll) 



represents the borderline separating the two regions in the definition of the function (|B10| ). 



It is clear from our derivation (and also easy to check) that / is continuous along this 
borderline. 



APPENDIX C: DISCUSSION OF THE FUNCTION / 

Since we do not have definite experimental values of c a and cp, but only bounds on 
these quantities (see Eq.( [3.17| )), which define allowed rectangles in the square < c a < 1, 
< cp < 1, we are interested in the behaviour of / in order to evaluate the unitarity bound. 

From the partial derivative of / in the region y > g(x), 

^ = ^k<x(2-2s-y), (CI) 
ox ox 

one can see that, at fixed y, the function / increases monotonously from x = to the point 
x = 1 — y/2, where the partial derivative in Eq.( |Cl| ) is zero. The points x = 1 — y/2 lie on 
the straight line y±(x) = 2 — 2x. In the range 1 — y/2 < x < 1 the function / decreases 
monotonously. Taking into account the symmetry f(x,y) = f(y,x), we see that at fixed x 
the function / increases monotonously from y = to the point y = l — x/2, where the partial 
derivative of / with respect to y is zero. These points lie on the straight line y2{x) = l — x/2. 
Beyond this line / decreases monotonously. Note that both straight lines lie in the range of 

/a- 
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Figure ^| shows a contour plot of the function f(x,y), together with the lines yi and y-i 
which intersect at the point 

x = y = \ . (C2) 

At this point both partial derivatives of / are equal to zero and therefore the point (C2) 
corresponds to the absolute maximum of /, given by 

/ max = / 2 (|0 =^^ 0.385. (C3) 

This number constitutes the absolute upper bound for \I a p\. 
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FIGURES 



FIG.|l[ Upper bound for the transition probability of LBL reactor u e 's into all possible 
states, 1 - for the SBL Am 2 in the range 10" 1 eV 2 < Am 2 < 10 3 eV 2 . The solid 

curve is obtained from the 90% CL exclusion plot of the Bugey v e — > v e experiment (see 
Eq.(f4.4|)). The dash-dotted and dash-dot-dotted vertical lines depict, respectively, the ex- 
pected sensitivities of the CHOOZ and Palo Verde LBL reactor neutrino experiments. The 
shadowed region corresponds to the range of Am 2 allowed at 90% CL by the results of the 
LSND experiment, taking into account the results of all the other SBL experiments (see 
Eq.Q). 

FIG.^. Upper bound for the probability of z/ M — > v e and z/ M — > v e transitions in LBL 
experiments. The solid curve is obtained only from the 90% CL exclusion plot of the Bugey 
v e — > v e experiment (see Eq.(|Q|)), whereas for the long-dashed curve the 90% CL exclusion 
plots of the BNL E734, BNL E776 and CCFR — > v e and — > v e experiments (see 
Eq.([4.6|)) and for the short-dashed curve the 90% CL exclusion plots of the CDHS f3"7fl 



and CCFR |38| — > experiments (see Eq. (|4.12|) ) have been used in addition. The 



dotted, dash-dotted and dash-dot-dotted vertical lines represent, respectively, the expected 
sensitivities of the LBL accelerator neutrino experiments KEK-SK, MINOS and ICARUS. 
The darkly shadowed region is allowed by the results of the LSND experiment, taking into 
account the results of all the other SBL experiments, in the case of LBL oscillations in 
vacuum. The two horizontal borderlines correspond to the limits ( |1.4j) for Am 2 and the 
left borderline corresponds to the lower bound in Eq.( |4.10| ). The lightly shadowed region is 
allowed by the results of the LSND and all the other SBL experiments in the case of LBL 
oscillations in matter. The solid curve constitutes also an upper bound for the probability 
of v P — > v T and v e — > v T transitions. 



FIG.[3[ Contour plot of the function f(x,y) given in Eq. (|5.2|) . The dotted line is the 
borderline g(x) = 2(1 — x)/(2 — x) between the regions where f = fi and / = The two 
solid lines represent the functions yi(x) = 2 — 2x and 2/2(2;) = 1 — x/2. 

FIG.f|. Upper bound for the parameter \I e/1 \ which characterizes the CP-odd asymmetry 
in the — > v e channel for the SBL parameter Am 2 in the range lO^eV 2 < Am 2 < 10 3 eV 2 . 
The solid curve represents the upper function in Eq.( |5.4p and is obtained from the 90% CL 
exclusion plot of the Bugey v e — > v e experiment. The dotted curve improves the solid curve 
where a° < a° e /2 (the lower function in Eq. (|5.4j )). It is obtained from the 90% CL exclusion 



plots of the Bugey v e — > v e experiment and the CDHS and CCFR — > and — > z/ M 
experiments. The dash-dotted curve is obtained from the 90% CL exclusion plots of the 
Bugey v e — > v e experiment and the BNL E734, BNL E776 and CCFR — > v e and z/ M — > v e 
experiments (see the upper function in Eq. (|5.8| )). The shadowed region corresponds to the 
range (|1.4|) of Am 2 allowed at 90% CL by the results of the LSND experiment. The solid 



curve represents also an upper bound for \I K 



CT 
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FIG.[5[ Upper bound for the parameter \I^ T \ which characterizes the CP-odd asymmetry 
in the — + v T channel. The solid curve is obtained from the 90% CL exclusion plots 
of the CDHS and CCFR — > and — > experiments (see Eq.(|5.5|)). The dotted 
curve is obtained from the 90% CL exclusion plots of the FNAL E531 and CCFR — > v T 
experiments (see Eq.( f5.9| )). The shadowed region corresponds to the range ( |1.4| ) of Am 2 
allowed at 90% CL by the results of the LSND experiment. 
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Figure |3| 
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